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The non-leptonic weak processes s + u — > u + d and u + d — > s + u are known to dominate the 
dissipation mechanism responsible for the viscosity of strange quark matter in its normal phase. The 
rates of such processes remain unknown for many color superconducting phases of quark matter. 
In this paper, we partially fill up the gap by calculating the difference of the rates of the two non- 
leptonic weak processes in four transverse spin-one color superconducting phases of quark matter 
(slightly) out of /3-equilibrium. The four phases studied are the color-spin locked phase, the polar 
phase, the planar phase and the yl-phase. In the limit of vanishing color superconducting gap, we 
reproduce the known results in the normal phase. In the general case, the rates are suppressed 
relative to the normal phase. The degree of the suppression is determined by the structure of the 
gap function in momentum space, which in turn is determined by the pairing pattern of quarks. 
At low temperatures, the rate is dominated by the ungapped modes. In this limit, the strongest 
suppression of the rate occurs in the color-spin-locked phase, and the weakest is in the polar phase 
and the yl-phase. 
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I. INTRODUCTION 

The interior of neutron stars is made of very dense 
baryonic matter. Currently our knowledge regarding the 
actual state of such matter is incomplete. One commonly 
accepted hypotheses is that the densest regions inside 
neutron stars are made of quark matter [![. Moreover, 
such quark matter may be a color superconductor 0, H[ . 
(For reviews on color superconductivity see for example 
Refs. H i, i, 0, H ©, M, EL El.) From the viewpoint 
of basic research, it is of fundamental importance to test 
this hypothesis empirically. 

The way to test the idea regarding the presence of 
quark matter inside stars is to make predictions regard- 
ing physics processes that affect observable features of 
stars and then test them against the stellar data. One 
class of physics properties that are substantially modified 
by the presence of color superconducting quark matter is 
related to the rates of weak processes. Such processes, for 
example, affect the cooling rates [lo| and the suppression 
of the rotational (r-mode) instabilities [HI in stars. The 
latter in particular is determined by the viscous proper- 
ties of dense matter [r| . 

Theoretically, the ground state of baryonic matter at 
very high density corresponds to the color-flavor-locked 
(CFL) phase of quark matter [l6[. In this phase, quarks 
of all three colors and all three flavors participate in spin- 
zero Cooper pairing on equal footing. The rates of the 
weak processes and some of their effects on the physical 
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properties of the CFL phase of quark matter have been 
discussed in Refs. 0, 111 EI SMI . 

With decreasing the density, the CFL phase should 
break up. This is due to the disruptive effects of a 
large difference between the masses of the strange quark 
and the light (up and down) quarks [2q . Such a differ- 
ence leads to a mismatch between the Fermi momenta 
of quarks and, therefore, spoils the "democratic" pair- 
ing of the CFL phase. When the CFL phase breaks up, 
another type of spin-zero color superconductivity, the so- 
called two-flavor color superconducting (2SC) phase [H[ , 
can still be possible. In the 2SC phase, strange quarks do 
not participate in pairing. Also, up and down quarks of 
one color remain unpaired. Some weak processes in the 
2SC phase and their effects on the physical properties 
have been studied in Refs. [13, Hll- 

It is important to mention that matter inside stars 
is neutral (at least on average) and in /3-equilibrium. 
Enforcing these two conditions affects the pairing be- 
tween quarks and may disrupt the usual formation of 
cross- flavor spin-zero Cooper pairs [22| . In this case, the 
ground state of matter can be in other forms, for example, 
such as stable variants of cr yst alline [26j , gapless [13, HH , 
or other exotic phases [29L |30T|. When spin-zero pairing 
cannot occur, the ground state can be in one of the spin- 
one color superconducting phases, in which same flavor 
quarks combine to form Cooper pairs [HI, [H, HI, H3, HH • 

Compared to the spin-zero case, the energy gap in spin- 
one color superconductors is likely to be about two orders 
of magnitude smaller. This means that the actual value 
of the gap may be somewhere in the range from 0.01 MeV 
to 1 McV. It appears that even such relatively small 
gaps can substantially affect the cooling rate of a quark 
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FIG. 1: Feynman diagrams of the non-leptonic weak processes 

star [36| . By the same token, such gaps can strongly 
modify the rates of the non-leptonic weak processes and, 
thus, affect the viscosity of stellar quark matter. 

The bulk viscosity in the normal phase of three-flavor 
quark matter is usuall y d ominated by the non-leptonic 
weak processes [IE HEM 0, SI El- (The correspond- 
ing processes are diagrammatically shown in Fig. [1]) It 
was argued in Ref. [43[, however, that the interplay be- 
tween the Urea and non-leptonic processes may be rather 
involved even in the normal phase of quark matter. In- 
deed, because of the resonance-like dynamics responsible 
for the bulk viscosity and because of a subtle interference 
between the two types of the weak processes, a larger 
rate of the non-leptonic processes may not automatically 
mean its dominant role. In fact, it was shown that the 
contributions of the two types of weak processes are not 
separable and that, at low frequencies relevant for some 
pulsars, taking into account the Urea processes may sub- 
stantially modify the result [43 |. 

Currently it remains unknown how a similar interplay 
between the two types of weak processes is realized in 
spin-one color superconducting phases. Primarily, this 
is because the rates of the non-leptonic weak processes 
in the corresponding phases have not been calculated. 
(Note that the rates of the Urea processes in several 
spin-one color superconducting phases were obtained in 
Refs. [IB, The purpose of this paper is to study the 
corresponding non-leptonic rates. 

The rest of the paper is organized as follows. The 
derivation of a general expression for the non-leptonic 
rate, based on the Kadanoff-Baym formalism [46J , is pre- 
sented in the next section. The structure of the quark 
propagators in spin-one color superconducting phases is 
described in Subscc. Ill Al This is used in Subsec. Hi Bl to 
derive the imaginary part of the lU-boson polarization 
tensor, which is the key ingredient in the expression for 
the rate. The net rate of the d-quark production (i.e., the 
difference of the rates of s+u — > u+d and u+d — > s+u) in 
the case of a small deviation from chemical equilibrium is 
obtained in Sec. lIIII There we also present the numerical 
results for each of the following spin-one color supercon- 
ducting phases: the CSL phase f Subsec. IIII A[) . the polar 
phase (Subsec. IIII Bp . the A-phase (Subsec. IIII C[) . and 
the planar phase fSubsec. IIII Dp . In Sec. IIVI we discuss 
the main results and their physical meaning. Two Ap- 
pendices at the end of the paper contain some details, 



used in the derivation of the rate. 
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FIG. 2: Feynman diagram for the d-quark self-energy. The 
particle four-momenta are shown in parenthesis next to the 
particle names. 



II. FORMALISM 

In order to calculate the rates of the non-leptonic pro- 
cesses, we use the same approach as in Refs. [HI, H(| 
IH, H, EH- It is based on the Kadanoff-Baym formal- 
ism (46[ . The starting point of the analysis is the general 
Kadanoff-Baym equation for the Green functions (prop- 
agators) of the down (or strange) quarks. After applying 
the conventional gradient expansion close to equilibrium, 
we derive the following kinetic equation for the d-quark 
Green function: 

ij^[ 70 S£(i\)] = -Tr[5>(P 1 )S<(P 1 ) 

-E>(P 1 )5<(P 1 )]- (1) 

Here we denote the quark four-momenta by capital let- 
ters, e.g., P = (po, p), where po is the energy and p is the 
three-momentum. The structure of the quark Green's 
functions S f< (P 1 ) and S > (Pi) in spin-one color super- 
conducting phases will be discussed in the next subsec- 
tion. To leading order, the quark self-energies E<(Pi) 
and E>(Pi) are given by the Feynman diagram in Fig. [2] 
This translates into the following explicit expression: 

S<,> (^) = ^/ 0r;_5<>(p 4 )«<(o), 

(2) 

where, by definition, M\y and Q = Pi — Pa are the 
mass and the four-momentum of the IF-boson, respec- 
tively. (Note that the large hierarchy between the W- 
boson mass and a typical momentum transfer Q < 1 MeV 
justifies the approximation in which the lU-boson prop- 
agator is replaced by 1/M^.) As seen from the diagram 
in Fig. [2l the expression for the polarization tensor of the 
IU-boson is given by 
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n^>(Q) 



(27T) 



-Tr 



C i+ 5>'<(P 2 + Q)r^_s<'>(p 2 ) 



(3) 



In the Nambu-Gorkov notation used here, the explicit form of the (tree-level) vertices for the weak processes d 
u + W~ and s <-> u + W~ reads |25| 



ud/us,i: 



e V ud/us ( ^(l-J 5 )T ud/uS: ± 

2V2sin6» w V 







(4) 



These are given in terms of the elements of the Cabibbo-Kobayashi-Maskawa matrix V uc i and V USl and the weak 
mixing angle 9w- By construction, the r-matrices operate in flavor space (u,d,s) and have the following form: 



T u d,+ 




Tud,- 




r us,+ 





By making use of Eqs. (JTJ) and ©, the kinetic equation takes the following form: 



dt 



dP 4 



tv[ 7 o^(Pi)] = -jjr J j^^[s^(Pi)Kd,-S^(P^Kd,+^u(Q) - <d,-S>(Pi)Kd, + sZ(Pi)iiUQ)\ ■ ( 6 ) 



(5) 



The physical meaning of the expression on the left hand 
side of this equation is the time derivative of the d-quark 
distribution function. By integrating this over the com- 
plete phase space, we obtain the net rate of the d-quark 
production: 



id f dPi 



idtj (2?r)' 



(7) 



Then, by making use of the kinetic equation ([T]), we de- 
rive 



4M^ 



d 4 Pi 



dP 4 



(2tt) 4 J (2tt) 4 

Tr {snPi)Kd,S<(PdKd, + ) n>,(Q) 
-Tr (r^_5>(p 4 )r^ + 5<(Pi)) n<,(Q)] . (8) 



This rate should be non-vanishing only if the rates of 
the two non-leptonic weak processes u + s — ► d + u and 
d + u — > it + s differ. In /3-cquilibrium, in particular, 
the latter two should be equal and the net rate of the 
d-quark production should vanish. The corresponding 
state of equilibrium in dense quark matter is reached 
when the chemical potentials of all three quark flavors 
are equal, i.e., fi u = Hd = Ms- (For simplicity, here it is 
assumed that all three quark flavors are approximately 
massless and, therefore, that the electrical neutrality of 
quark matter is achieved without the need for the elec- 
trons.) 

When the system is forced out of equilibrium, e.g., dur- 
ing the density oscillations caused by the collective modes 
of stellar matter, small deviations from /3-cquilibrium are 
induced. For our purposes, the corresponding state can 



be described by the following set of the chemical poten- 
tials: n u — fid — fi- and fi s = fi+ Sfi, where 5[i is a small 
parameter that characterizes the magnitude of the de- 
parture from the equilibrium state. Out of equilibrium, 
the net production of d-quarks may be nonzero. For ex- 
ample, if Sfi > (Sfi < 0) the system has a deficit (an 
excess) of the down quarks and an excess (a deficit) of 
the strange quarks. Then, one of the weak processes, i.e., 
u + s — > d + u (d + u — > u + s), will start to dominate over 
the other in order to restore the equilibrium. The net 
rate of the d-quark (or equivalently s-quark) production 
characterizes how quickly this happens. 

In order to calculate the net rate of d-quark produc- 
tion, however, one needs to know the explicit structure of 
the quark propagators in the specific spin-one color su- 
perconducting phases. The knowledge of the quark prop- 
agators is also needed for the calculation of the polariza- 
tion tensors H<„(Q) and U^(Q). These are discussed in 
the next two subsections. 



A. Quark propagator in a spin-one color 
superconductor 

In general, Cooper pairs in spin-one color supercon- 
ducting phases are given by diquarks in a color antitriplct 
(antisymmetric) and spin triplet (symmetric) state. De- 
pending on a specific color-spin structure, which is deter- 
mined by the alignments of the antitriplet in color space 
and the triplet in spin (coordinate) space, many inequiv- 
alent color superconducting phases may form. 

Each phase is unambiguously specified by the struc- 
ture of the gap matrix, which is commonly written in the 
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following form 



*(P)" 



e=± 



(9) 



where 4> e (P) is the gap function. The Dirac matrices 
Ap = (1 + e7o7 ■ p)/2, with e = ±, arc the projectors 
onto the positive and negative energy states. The color 
structure of $(P) + is determined by 



J 4 A, 



p J COS i 



+ 7^ sin 6 



(10) 



where (J i y k = —ie 1 ^ are the antisymmetric matrices in 
color space, p = p/p is the unit vector in the direction 
of the quasiparticle three-momentum p, and 7^ = "f J — 
P'il ' p)- The explicit form of the 3x3 matrix Ay and 
the value of the angular parameter 9 determine specific 
phases of superconducting matter. Among them, there 
is a number of inert and noninert spin-one phases [35| . 
which arc naturally characterized by the continuous and 
discrete symmetries preserved in the ground state. 

In the two special cases, — and 8 = ir/2, the corre- 
sponding phases are called longitudinal and transverse, 
respectively. In this paper we focus on the transverse 
phases (6 = 7r/2), in which only quarks of the opposite 
chiralities pair and which have lower free energies than 
the longitudinal phases [32] ]. To further constrain the 
large number of possibilities, we concentrate only on the 
following four most popular ones: the color-spin locked 
phase (CSL), the A-phase, the polar phase and the planar 
phase. 



The structure of the matrices A,- 



and A4 P for the 



mentioned four phases are quoted in the first two 
rows of Tab. Q] (for more details see Refs. [HI). In 
the corresponding ground states, the original symme- 
try SU(3) e xSO (3) j xU(l)em of one-flavor quark matter 
breaks down to M H \M M, HI 



SO(3)j _ 
SU(2) C x SO(2)j x 17(1) 
SU(2) C x SO(2)j x U(l) 
SO(2)j x £7(1) cm 



cm 
cm 



(CSL), 
(A-phasc), 
(polar), 
(planar), 



respectively. 

In spin-one color superconductors, there is no cross- 
flavor pairing and, therefore, the quark propagator is di- 
agonal in flavor space, i.e., 



S(P) = diag[S u (P), S d (P), S S (P)}. 



(11) 



The Nambu-Gorkov structure of each flavor-diagonal el- 
ement is given by 



Sf>>(P) 



G<>>(P) F<>>(P) 
F<>>(P) G<>>(P) 



(12) 



where f = u,d, s. The normal (diagonal) and anomalous 
(off-diagonal) components of the Nambu-Gorkov propa- 
gator have the following structure [32| : 



G<;>(P) = 1oA*52v± r G<>> f (P).. 



(13) 



Ffi?(P) = -7oA4p7oE 7 'p^ A p e ^( P )' (I 4 ) 



(15) 



Here, r labels different quasiparticle excitations in color- 
superconducting quark matter. The matrices V~ r and 
V pr are the projectors onto the subspaces spanned by 
the eigenvectors of M P M P and j a M p M p j° , respec- 
tively. The explicit form of the projectors for each phase 
can be found in Ref. [36[. It should be noted that both 
matrices M P M P and 7°A4pA / J P 7 have the same set of 
eigenvalues A PjI ., 

M P M P = EVrPp.r , (16) 

r 

7°X P X P 7° = ' 

r 

The list of all eigenvalues as well as their degeneracies 
are given in the last three rows of Tab. [J] Each of the 
eigenvalues determines a quark quasiparticle with the fol- 
lowing dispersion relation: 



c p,r 



/ = \/(.P-H-) 2 + WK,rJ- (18) 



The separate components of the propagators in subspaces 
spanned by the eigenvectors, see Eqs. (fl3|) . (fT4|) and (|15jl . 
can be conveniently rewritten in terms of the correspond- 
ing distribution functions f(e P ,r,f) and the Bogoliubov 
coefficients B pr f, i-e., 

G±,r,f{P) = -2tti B p,r,//( ee P,r,/)<5(po ± [if - eep, r , f ), 

e=± 

(19) 

G±, r ,f(P) = -2tti ^2 £p,r,//(-«ip,r,/)<S(P0 ± M/ ~ ee P,r,/), 

(20) 

F ±,r,f(P) = 27Ti ~ e f( e€ P*rj) S (P0 T l*f ~ eCp.r,/), 

Ze P,r,f e __j_ 

(21) 

Ze P,r,f e _ ± 

(22) 

The Bogoliubov coefficients and the fermion distribution 
function are defined as follows: 



B 



P,r,f 
/(e) 



1 

cxp(y) + 1 ' 



(23) 
(24) 



5 





CSL phase 


planar phase 


polar phase 


yl-phase 


Ay 


5y 


SilSjl + 8i2&j2 




Sis(8ji + i S j2 ) 


M p 


J ' 7± (p) 


Ji7-L,i(p) + ^27±,2(P) 


J r 37±,3(p) 


J3[tx,i(p) + *7-L,2(p)] 


A p ,i (m) 


2(8) 


1 + cos 2 <? p (8) 


sin 2 P (8) 


(l + |cos0 p |) 2 (4) 


A p ,2 (n 2 ) 


0(4) 


0(4) 


0(4) 


(l-|cose p |) 2 (4) 


A P ,3 (n 3 ) 








0(4) 



TABLE I: Matrices Ay and M p as well as the eigenvalues A p , r with the corresponding degeneracies n r in four spin-one color 
superconducting phases. The angle between p and the z-axis is denoted by P . 



[Note that /(— e) = 1 — /(e).] The quark propagators in 
Eq. (JT3J) can now be used to derive the general expres- 
sions for H< U (Q) and II^, (Q). This is done in the next 
subsection. The results are then used to calculate the 
rate T d in Eq. ©. 



B. W-boson polarization tensor 

The VF-boson polarization tensor is given in terms of 
the quark propagators in Eq. By taking into account 
the Nambu-Gorkov and flavor structure of the weak inter- 
action vertices in Eq. ([4]), as well as the quark propagator 
in Eq. (jT2J), we derive 



ie 2 V 2 

n<j>(Q) = - 



xTr 



8sm z 6 w J (2tt) 4 
7*(1 - 7 5 )G> + < (P 3 )7^(1 - 7 5 )G<;>(F 2 ) 



+7^(l+7 5 )G^(P 3 )7 ,y (l+7 5 )G s < ;L > (P 2 )j , (25) 

where we introduced the notation P 3 = P 2 +Q. Note that 
the anomalous (off-diagonal) elements of the Nambu- 



Gorkov propagators dropped out from the result. This 
is the consequence of the electric charge conservation. 
In calculations, this comes about as a result of the spe- 
cific flavor structure of the weak interaction vertices in 
Eq. ([4]). One can further simplify the result for the polar- 
ization tensor in Eq. (|25[) by noticing that the two terms 
on the right hand side are equal. From physical view- 
point, this is related to the fact that the two terms are 
the charge-conjugate contributions of each other. After 
taking this into, we arrive at the following expression for 
the polarization tensor: 



n^>(Q) 



4 sin 2 e 



(2tt) 



-Tr 



7 M (l-7 5 ) 



xG> : <(P 3 )y"(l-y 5 )G<;>(P 2 ) . (26) 



Then, by using the explicit structure of the normal com- 
ponents of the u- and s-quark propagators, defined in 
Eq. flU]) , we obtain 



n<->(Q) 



4 sin 2 (9 



w 



2tt) 



rTr 



'(1 - 7 5 )7oA p3 £>+ ^X^AWi 1 - 7 5 )7oA p2 , r2 G<;> iU (P 2 ) 



This can be rewritten in an equivalent form as 



where, by definition, the tensor (p,p') is given by the following trace (in color and Dirac spaces): 

r;;'(p,p') = Tr[ 7 ^(l - 7 5 hoA p P+ r7 l '(l - 7 5 )7oA p ,P+ r ,\ . 



(27) 
(28) 



(29) 



This trace was calculated for each of the four spin-one color superconducting phases in Ref. [361 ]. For convenience, 
the corresponding results are also quoted in Appendix [Aj 

Finally, by making use of Eqs. (TIT))) and PO)) . we arrive at the following expression for the VF-boson polarization 



6 



tensor: 



n <,> 



i7Te 2 V„ 



us 



2 sin 2 9 



d :i 



P2 



W - 

r2,r 3 ,ei,e2 

x /(±eie P 3,r 3 ,s)/(Te2ep 2 ,r 2 , tl )<5('7o + 8[i - e x e P3 



(2tt) 3 



e 2 

P2,'"2,M 



r< m>^ 



(30) 



Here we denote 8 [i = \i s — \i u and assume that the upper (lower) sign corresponds to Il < (n > ). It should be mentioned 
that one of the 8- functions was used to perform the integration over p2,o- 



III. CALCULATION OF THE RATE 



In this section we derive a general expression for the net rate of the d-quark production in spin-one color super- 
conducting quark matter close to chemical equilibrium. The corresponding rate is formally defined by Eq. ([8]). By 
making use of the quark propagators and the IF-boson polarization tensor, derived in the previous section, we obtain 



16M^ sin 2 9 W 



f d 4 P 1 f d 4 P 4 ^ 
J (2^7 (2^^ T - l(P4 ' Pl 



1-1,1-4 
> 



where we used the following results for the traces: 



Tr(5>(P 1 )r^_5<(P 4 )r^, + 
Tr(r^_5>(P 4 )0, + 5<(P 1 ) 



e 2 V 2 



4 sin 2 9 W 



e 2 V 



4 sin" 2 6 W 



ff; E Tf"r 1 {i > iii > i)G%, riyU {Pi)G+ rud (Pi). 



(31) 

(32) 
(33) 



As in the calculation of the polarization tensor, the anomalous (off-diagonal) Nambu-Gorkov components of quark 
propagators did not contribute to these traces. This is the consequence of the specific flavor structure of the weak 
interaction vertices . 

After making use of Eqs. (JTSJ), ([20]) and ((201), we obtain 



T d = 2 7 n 4 G 2 F V 2 d V 2 s 



E E 



d 3 px d 3 p 2 d 3 p 3 d 3 p 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 



(1 - Pl • P2)(l - P3 ' p4)Wr 1 r 1 (p4,PlVr 3 r 2 (p3,p2) 



X B p 1 1 ,r 1 ,d- B P2,r2,« B P3,r3, S B p 4 4,r 4 ,u'5(Pl + P2 P3 - Pi)8(ei £ pi , ri , d + e 2 e p2 ,r 2 ,u ~ e 3 e P 3,r 3 ,s ~ e 4 <W 4 ,« + Sfl) 

x [/(eie piiI . l!( z)/(e2ep2,r2,ti)/(- e 3ep3,r3,s)/(- e 4ep4,r4,«) ~ /(— eie pi ,ri,d)/(— e 2%>2,r 2 ,i0/( e 3 e P3,i-3,s)/( e 4ep4,i-4,«)] • 

(34) 



r 



In derivation, we used the definition of the Fermi constant 
in terms of the VT-boson mass, 

e 2 

G F = = (35) 

4^2 sin 2 9 W M^ 

and the following Lorentz contraction: 

^(p^,pd)?;r 2 (p3,p2) = 16(1 -pi • P2 ) 

x(l - p 3 ■ p4)w r4ri (p4,Pl)w r3 r 2 (P3,P2) ,(36) 

where w rr <(p,p') denotes a color trace that involves a 
pair of quasiparticles (r and r') with the given directions 
of their three-momenta (p and p') in a specific spin-one 
color superconducting phase. The corresponding traces 
for all four phases are listed in Appendix [Al 



Formally, the expression in Eq. (|34|) gives the net rate 
of the d-quark production in quark matter away from 
chemical equilibrium. The first term in the brackets de- 
scribes the production of (i-quarks due to s + u — > u + d, 
while the second one describes the annihilation of d- 
quarks due to u + d ^ s + u. 

Here it might be instructive to note that the above 
expression for the rate T d resembles the general result 
for the net rate of the <i-quark production in the normal 
phase of strange quark matter [47|. The key difference 
comes from the presence of the Bogoliubov coefficients 
Bp.r.f and the u> rr i(p, p') functions that account for a 
non-trivial quark structure of the quasiparticles in spin- 
one color superconductors. Naturally, when such quasi- 
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particles are the asymptotic states for the weak processes, 
the amplitude is not the same as in the normal phase. 

The degree of departure from /3-equilibrium and, thus, 
the net rate is controlled by the parameter 5fj, = fi s — \i d - 
When 5[i = 0, the expression in the square brackets of 
Eq. (l34|) vanishes and T d = 0. When Sfj, ^ 0, on the 
other hand, one has 

T d ~ \6fi (37) 



to leading order in small S/j, J48|. Note that the overall 
sign was chosen so that A is positive definite. (Recall 
that a positive 5fi means an excess of strange quarks, 
which should drive a net production of d-quarks, while a 
negative S/i means a deficit of strange quarks, which will 
be produced by annihilating some d-quarks.) 

From the general expression in Eq. (|34| . we derive 



5A 

A — 



E E J d 3 p 1 d 3 p 2 d 3 p 3 d 3 p 4 (l - pi ■P2)(l-P3-P4)w r4ri (p4,pi)w r3r2 (p3,p2) 



2 11 n 5 u 5 T 

rir2r3r 4 eie 2 e 3 e 4 

x • B pi,n,d B P2^2,«- B p 3 3,''3,s- B p 4 4^4,u (5 (Pi + P2 ~ P3 - P4)£(ei e Pl >ri , d + e 2 e P2 ,r 2 ,« - e 3 ep3 !r3 , s - e 4 e P4ir4iU ) 

x /(-eiep ljri>d )/(-e2ep 2>r2) „)/(e 3 ep 3ir 3 jS )/(e4ep 4ir4j „). (38) 



where 



A = F V K "V 5 T 2 (39) 
is the corresponding A-rate in the normal phase of strange quark matter |4?| . 



A. Analysis of the rate in CSL phase 



Out of the four spin-one color superconducting phases studied in this paper, the CSL phase is special. This is the 
only phase in which the dispersion relations of quasiparticles are isotropic. As a result, the corresponding rate is the 
easiest to calculate. In this subsection, we analyze the A-rate in the CSL phase in detail. 

Let us start by noting that the explicit form of the w rr '(p, p')-functions in the CSL phase is given by 

wii(p.p') = l + i(l + p- P ') 2 , (40) 

wi 2 (p, p') = «ai(p, p') = 1 - i(l + p • P ') 2 , (41) 

^2 2 (P,P') = ^(1+P-P') 2 - (42) 

(See Appendix [Al and Ref. [3^|.) By making use of these expressions and introducing the following notation for the 
angular integrals: 

Fr ir2 r 3 r4 = J dQl Jdil 2 J dfl 3 J dfli (1 - f> 3 ' Pi)(l ~ Pi • P2)w r4 n (P4, Pl)^r 3 r 2 (p3» P2)5(Pl + P2 ~ P3 ~ P4) , (43) 

we arrive at the following representation for the A-rate in the CSL phase: 
5A /i 3 



A (CSL) 



pOC poo poo poo 

Jo Jo Jo Jo 



x /(-eie P1 ,n,d)/( )/(e3e P 3,r 3 ,s)/(e4ep 4 ,r4,uM( e i e pi,ri,d + e 2£ P2 ).(44) 

I 



Here we took into account that, to leading order in in- 
verse powers of fj,, the absolute values of the quark three- 
momenta can be approximated by fi. In the same approx- 
imation, the explicit form of functions F rir2rari are given 



in Appendix |B"1 All of them are proportional to 1/n 3 . 
This factor cancels out with the overall fj, 3 in Eq. (j44|) . In 
order to perform the remaining numerical integrations, it 
is convenient to introduce new dimcnsionlcss integration 
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variables Xi = {pi — [t)jT instead of pi (i = 1,2,3,4). The 
integration over 24 is done explicitly by making use of the 
(5-function. The remaining three-dimensional integration 
is done numerically, using a Monte-Carlo method. One 
finds that the ratio A ( - CSL - ) /Ao is a function of a single 
dimensionless ratio, 4>/T. 

Before proceeding to the numerical results, it is in- 
structive to analyze the limiting case of low temperatures 
(or alternatively very large 4>/T). In this limit, only the 
ungapped r = 2 quasiparticlc modes should contribute 
to the rate. The corresponding contribution is easy to 
obtain analytically, i.e., 



A (CSL) 



7? 

2i- 3 r 4 r r 1 r 2 r 3 r i 



-^-A w 0.034A , for - -> 00. (45) 
27027 T \ ) 



The subleading correction to this result is suppressed by 
an exponentially small factor exp(— y2(j)/T). (Note that 
y/2 in the exponent is connected with the conventional 
choice of the CSL gap, which is V20 rather than </>.) 

It might be instructive to mention that the asymptotic 
value in Eq. (|4"5|) is substantially smaller than Ao/9, which 
is the corresponding contribution of a single ungapped 
mode in the normal phase. The additional suppression 
comes from the functions W22(p4,pi) and W22(p3,p2) 
which modify the amplitude of the weak processes with 
respect to the normal phase. Except for the special case 
of collinear processes (i.e., p4 parallel to pi and p3 paral- 
lel to P2), the corresponding w-functions are less than 1, 
see Eq. (|4"2"j) . Interestingly, this kind of suppression is a 
unique property of the non-leptonic rates and is not seen 
in analogous Urea rates because the latter are dominated 
by the collinear processes [H, HH . 

All our numerical results for the A-rates as a function 
of <j)/T are shown in Fig. [51 [49jj . In the case of the CSL 
phase (black points and the interpolating line in Fig. [3]), 
we used the Mathematica's adaptive quasi-Monte-Carlo 
method to calculate the A-rate. In order to improve the 
efficiency of the method, we partitioned the range of in- 
tegration for each of the three dimensionless integration 
variables Xi = (pi — fJ.)/T into several (up to 6) non- 
overlapping regions. This approach insures that the main 
contribution, coming from a close neighborhood of the 



Fermi sphere, is not lost in the integration over a for- 
mally very large phase space. 

As seen from Fig. [3l the numerical results smoothly 
interpolate between the value of the rate in the normal 
phase Ao and the asymptotic value of the rate due to the 
CSL ungapped modes, given by Eq. (|4"5|) . 



B. Analysis of the rate in polar phase 

Unlike the CSL phase, the polar phase is not isotropic. 
However, it is the simplest one among the other three 
phases. While the dispersions relations of its quasipar- 
ticles de pend on the angle 9 p between the momentum p 




0/T 



FIG. 3: (Color online) Numerical results for the A-rate in four 
different phases of spin-one color superconducting strange 
quark matter [4S|]. The error bars show the statistical error 
estimates in the Monte-Carlo calculation of the rates. The 
horizontal dashed lines correspond to the contributions of the 
ungapped modes in the limit of large <j>/T (or equivalently the 
limit of low temperatures) . 



and a fixed z-drection, its u> rr >(p, p')-functions are inde- 
pendent of the quasiparticle momenta, i.e., 



U) rr '(p,p') = n r 8rr> , 



(46) 



with n\ = 2 and 712 = 1, sec Appendix lAl Taking this 
into account, the corresponding A-rate takes a simple 
form: 



5 A f 

A (P ° lar) = — e Y. n n n r2 / rf 3 Plrf 3 P2d 3 P 3 rf 3 P4(l-pl-p 2 )(l-p3-p4)'5(Pl+P2-P3-P4) 



2 n 7rVT 3 

I J_ I ^ l-l 1 -- i t 

X 5 p 1 ,r 1 ,d S P2,r2 1 « 5 plr2,s 5 pl,ri,«/( _e l e Pi,ri,d)/(-e2ep 2 ,r 2 ,«)/(e3ep 3 ,r 2 , s )/(e4e P4 ,r 1 ,«^ 
X ^(ci6p i ri d -j- C2^ P 2,r 2 ,ii ^3^p3,r2,s ^4^p4,ri,u)- 

By making use of the first 5-function, we easily perform the integration over P4. We can also perform the integration 
over one of the remaining polar coordinates. This is possible because the integrand depends only on the two indepen- 
dent combinations of the polar angles, i.e., (f>\ = <p\ — ^3 and y5 2 = ^2 ~ <P3- By using (p\ and (pi as new integration 
variables (for simplicity of notation, the tildes are dropped in the following) , we see that the integrand is independent 
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of the variable 953. Finally, by approximating p\p\p\ — /i 6 in the integration measure, we rewrite the expression for 
the rate as follows: 



2 18 7r 4 T 

nr 2 



1 



2tt r 2ir 



r\ roc poo />oo />! />! />! 

A(P° lar > ~ -^i-^n ri n r2 / dx x / dx 2 / dx 3 / d& / <*6 / #3 / ^1 / #2 



(1-C0S6» 12 )(1-C0S6»3 4 ) Dei n e 2 D e 3 n e4 



Dei TDe2 dC3 ne 

co S h(ie xliSl)ri )co S h(Ie X2ie2ir2 )cosh(ie X 3 i€3)r2 )co S h(ie X4ie4)ri ) eie ^ e4 x ^ *^ 

x 8(eie Xl £ 1>ri + e2£ X2 ,Z2,r2 ~ e 3 e x 3 ,i 3 ,r 2 - e4e X4i g 4jri ), (48) 



where the new integration variables are Xi = (pi — /i) /T 
and 6 = cos 6 Pi . By definition, the dimensionless energy 
is 



(49) 



with Af,x = 1 — 6 an d A^ 2 = 0, and the new Bogoliubov 
coefficients arc 



1 t x 



(50) 



Note that the expressions for X4, 6 and cos #34 in Eq. (|48|) 
are given by 



X4 = 



Pi - A' 
T ' 

+ P26 -Ks6 



P4 

pi COS 6<i3 + p 2 COS 6*23 - P3 
COS0 34 = , (51) 

Pi 

where P4 = |pi + P2 — P3I is a function of Xi (i = 1, 2, 3) 
and the three cosine functions, 



cosflia = 66 + V 1 _ ft V 1 ~ ?2 cos(^x - ip 2 ) 



cos 6»i3 = 66 + yi-C 2 yi-C3cos(^i), 

cos023 = 66 + V 1 " ^2^1-^1 cos( V2 ). (52) 

In the calculation, we used a customized Monte-Carlo 
method in order to improve the statistical error of the 
integration over Xi (with £ = 1,2, 3). To this end, we used 
a special type of importance sampling, which is motivated 
by the fact that the main contribution to the rate should 
come from the region near the Fermi surface. In order to 
implement this, we utilized random variables distributed 
according to the Gaussian distribution [50| : 



P(x) 



V2 



1 / (x-Xq) 2 

exp 1 



7T(7^ 



2a 2 



(53) 



where xq and a are the mean and the width of the dis- 
tribution, respectively. This was applied to the numer- 
ical integration over the dimensionless variables a;, = 
(pi — \i)jT (i = 1,2,3), in which case we took xq = 



and cr = 3. In order to generate independent variables 
(e.g., x\ and X2), distributed according to Eq. (|53"|) , we 
applied the Box-Mullcr transform, 



x\ = x + ay/— 2 \nu\ cos(27TW2), (54) 

X2 = %o + en/— 21nuisin(27TU2), (55) 

where u\ and 1*2 are two independent variables, uniformly 
distributed in the range from to 1. 

In our numerical calculation, we also used a Gaussian 
function to approximate the (5-function responsible for 
the energy conservation in the expression for the rate 
(j48]l . For this purpose, we used the width of the distri- 
bution do = 0.2. This appeared to be sufficiently small to 
avoid strong violations of the energy conservation in the 
weak processes and, at the same time, sufficiently large 
to use in a Monte-Carlo integration with the number of 
(eight-dimensional) random points on the order of 10 6 (in 
a Mathematica code) or 10 7 (in a Fortran/ C++ code). 

The numerical results for the A-rate in the polar phase 
are shown by the blue squares (and the interpolating 
line) in Fig. [3] At vanishing 4>/T, the rate coincides with 
that in the normal phase. At asymptotically large value 
of <j)/T, on the other hand, the rate approaches Ao/9. 
This value is marked by the purple dashed line in the 
figure. Theoretically, the rate is dominated by the un- 
gapped modes (n = r 2 = 2) in the <f>/T — > 00 limit. The 
corresponding contribution can be obtained by analytical 
methods as follows. We start by pointing that the Bo- 
goliubov coefficients for the ungapped modes are equal 
to the unit step functions: B^, ^ 2 = @( — e i x i)i where by 
definition Q(x) = 1 for x > and Q(x) = otherwise. 
Since these Bogoliubov coefficients are nonzero only for 
ej = sign(— Xi), each sums over ej effectively reduces to 
a single contribution. By taking this into account and 
making use of the result for the angular integration, Kq, 
defined in Appendix [Bl we derive 

\ />oo />oo poo poo 

A inSd ~ fiZS dx ^ dx * / dx 3 / 

J — 00 J — 00 J — 00 J — 00 

S(— x\ — x 2 + x 3 + X4) 



(e^ 1 + l)(e x i + l)(e~ x 3 + l)(e~ x * + 1) 



(56) 



It should be noted that the numerical results for the po- 
lar phase in Fig. [3] approach this asymptotic value very 
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slowly. We can speculate that this indicates a weak 
(probably, power-law) suppression of the contribution of 
the gapped (mixed with ungapped) modes to the rate. 
The key feature responsible for this behavior in the po- 
lar phase is the presence of gapless nodes at 6 p = and 
8 P = 7r in the dispersion relation of the gapped modes. 
As we shall see below, the same qualitative property is 
shared by the ^4-phase, whose gapped modes also have 
a node at 8 p = tt. In contrast, the rates in the CSL 
and planar phases, whose gapped modes have no gap- 
less nodes, show asymptotes that are consistent with the 
rapid, exponential approach to their asymptotic values. 

C. Analysis of the rate in yl-phase 

The analysis in the j4-phase of spin-one color supercon- 
ducting matter can be performed along the same lines as 
in the polar phase. The apparent complication of the 
yl-phase is the existence of three, rather than two dis- 
tinct quasiparticle excitations. However, it appears that 
the contributions of the two gapped modes (r = 1,2) 
can be replaced by a single contribution of a modified 
mode with the energy e p = yj(p — /i) 2 + |0| 2 A p where 
A p = (1 + cos6* p ) 2 (cf. the dispersion relations of the 
modes r = 1,2 in Tab. UJ. This alternative representa- 
tion is possible because of the special, separable structure 
of the corresponding u> rr i (p, p') functions in the ^4-phasc. 
As seen from the expressions in Eq. (|A7|) . the mode r = 1 
contributes only when cos 9 P of the corresponding quasi- 
particle is positive, while the mode r — 2 contributes 
only when cos# p is negative. Then, when the contribu- 
tions are nonvanishing, one always gets w rr /(p,p') = 2. 
By also noting that the corresponding eigenvalues 

A p ,i = (1 + |cos6»p|) 2 for cos(9p>0 (57) 

and 

Ap.2 = (1 - Icos^pl) 2 for cos6» p <0 (58) 

formally take the same form, i.e., A p = (1 + cos6> p ) 2 , we 
conclude that the sum over the original modes r — 1, 2 in 
the rate can indeed be replaced by a single contribution 
of the modified mode as defined above. By making use 
of this observation, the general expression for the rate in 
the ^4-phase takes the form, which is similar to that in the 
polar phase, see Eq. (|47|) . but with a different dispersion 
relation of the (modified) gapped mode. 

By using a Monte-Carlo algorithm as in the previous 
case, we perform a numerical calculation of the A-rate in 
the ^4-phase. The corresponding results are shown by red 
stars (and the interpolating line) in Fig. [3] In the limit 
of large 4>/T, the rate is saturated by the contribution of 
ungapped modes, which is the same as in the polar phase, 
namely Ao/9. The derivation of this asymptotic expres- 
sion is the same as in the polar phase. The corresponding 
value is marked by the purple dashed line in the figure. 
A slow (probably, power-law) approach of the asymptotic 



value at <p/T — > oo is again associated with the presence 
of a gapless node (at 9 p = tt) in the dispersion relation 
of the (modified) gapped quasiparticles. 

D. Analysis of the rate in planar phase 

The calculation of the rate in the planar case requires 
the largest amount of computer time. One of the main 
reasons for that is the much more complicated expres- 
sions for the ui rr i(p, p')-functions (see Appendix!^}. The 
numerical results for the A-rate in the planar phase are 
shown by green diamonds (and the interpolating line) in 
Fig. [3] The asymptotic value of the rate at large <fi/T 
was extracted only numerically. By taking into account 
possible systematic errors (e.g., due to the overall nor- 
malization of the rate that may differ by up to 15% from 
the analytical estimate in the normal phase), we es- 
timate A( planar ) ~ (0.038 ± 0.003)A for <j>/T -> oo. Note 
that this is smaller than Ao/9, which is the contribu- 
tion of a single mode in the normal phase. As in the 
CSL phase, in the planar phase the additional suppres- 
sion comes from the w-functions for the ungapped modes. 



IV. DISCUSSION 

In this paper we derived the near-equilibrium rates of 
the net d-quark production (or equivalently the A-rates) 
due to the non-leptonic weak processes (i.e., the differ- 
ence of the rates of u + d — > u + s and u + s — » u + d) in 
spin-one color-superconducting strange quark matter at 
high density. The main numerical results are presented 
in Fig. O 

In the limit of 4>/T = 0, which is same as the normal 
(unpaired) phase of strange quark matter, our results re- 
produce the known result of Ref. [47[ ■ The effect of color 
superconductivity is to suppress these rates. The degree 
of the suppression depends on the details of the specific 
spin-one phases. To large extent, this is controlled by the 
value of the energy gap (more precisely, <fr/T) as well as 
its functional dependence on the direction of the quasi- 
particle momentum. At very large <p/T (or equivalently 
in the limit of low temperatures), the A-rates approach 
fixed values, which arc determined by the contribution of 
the ungapped modes alone. The corresponding limiting 
value is the smallest in the CSL phase. It is less than a 
third of the "canonical" value Ao/9 due to a single un- 
gapped mode in the normal phase of matter. The ad- 
ditional suppression comes from the modification of the 
quasiparticles due to color superconductivity. A similar 
observation applies to the planar phase. The rates in the 
other two phases, i.e., the polar and ^4-phase, approach 
the asymptotic values equal to Ao/9. 

The numerical results for the A-rates in Fig.[3]also indi- 
cate that the asymptotic approach to the limiting values 
can be qualitatively different in spin-one color supercon- 
ducting phases. In the case of the polar and A-phase, 
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the approach seems to follow a power law. In contrast, 
the approach appears to be exponential in the case of the 
CSL and planar phase. This qualitative difference can be 
easily understood. The power law is the consequence of 
the presence of gapless nodes in the dispersion relations 
of the gapped quasiparticlcs in the polar and ^4-phasc 
(the nodes are located at 9 P = and 9 p = it in the polar 
phase, and at 6 p = it in the ^4-phase). In the CSL and 
planar phase, the approach to the asymptotic value at 
4>/T —> oo is exponential because no gapless nodes are 
found in their gapped quasiparticles. (Note that a sim- 
ilar observation regarding the rates of the semi-leptonic 
processes was made in Ref. [36l |44||.) 

The results for the rates of non-leptonic weak pro- 
cesses, presented here, is an important ingredient for 
the calculation of the bulk viscosity of spin-one color- 
superconducting strange quark matter. If such matter 
is present inside neutron stars, its viscosity will be one 
of the mechanisms responsible for damping of the stellar 
r-mode instabilities [14j . 
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In the polar phase, the w rr '(p, p') functions do not 
depend on the quark momenta. They are given by the 
following expressions: 

cJn(p,p') = 2, (A4a) 
w 22 (p, P ') = 1, (A4b) 
wi 2 (p,p') = w 2 i(p,p') = 0- (A4c) 

In the planar phase, the explicit form of the Lu rr i (p, p') 
functions reads 



wii(p,p') = ^ I 3 + ? ?(P>P')] 



(A5a) 



wia(p,P ; ) - w 2 i(p,p') = 2-[l-»7(P,p')]> ( A5b ) 



W22(P,P') = ^ I 1 + ? ?(P>P')] 



(A5c) 



where 



»?(P,P 



, _ ipzP'z + (PxP'x + PyPyY - {PxP'y - PyP'xf 



[i + (p z ni + (P> z y 



. (A6) 

In the ^4-phase, there are three different quasiparti- 
cle branches (r = 1,2,3). Consequently, there are more 
u> rr '(p, p') functions, i.e., 

wn(p,p') = i[l + sgn(^)][l + sgn(^)], (A7a) 
w 22 (p,p') = hi - Bgn(p,)][l - sgn(&)], (A7b) 



APPENDIX A: COLOR AND DIRAC TRACES 

In this appendix, we write down the explicit expres- 
sions for the tensor (p, p'), defined by Eq^ (|2"9]) . The 
corresponding results were obtained in Ref. [361 ] - In gen- 
eral, one finds that 



^(p, P ')=T^( P ,p'K r ,(p,p'), 



(Al) 



where w rr '(p, p') are the functions determined by a spe- 
cific color-spin structure of the gap matrix, and 

T""(p,p') = Tr D [y(l - 7 5 )7°A p7 ' y (l - 7 5 )7°A p , 

(A2) 

The explicit form of all the components of this tensor can 
be also found in Ref. [36| . It is more important for us 
here to note that the following result for the contraction 
of this tensor with itself is valid: 

T^(p4,Pi)V(p3, P2 ) = 16(1 - Pl • P2 )(l - P3 • p 4 ) . 

(A3) 

Since an essential information regarding spin-one color- 
superconducting phases is carried by the u) rr >(p, p') func- 
tions, we also quote them here. (For more details, see 
Ref. [11].) 



wia(p,p') = -[l + sgn(^)][l-sgn(^)], (A7c) 
w 2 i(p,p') = \[l -sgn(p 2 )][l + sgn(&)], (A7d) 



wi 3 (p,p') = w 3 i(p,p') = 0, 

W 2 3(P,p') = W3 2 (p,p') = 0, 
W 3 3(P,P') = 1. 



(A7e) 
(A7f) 
(A7g) 



Finally, in the CSL phase, the corresponding functions 
are 



«n(p,p') = l + i(l+p-p') 2 , 



(A8) 



«i 2 (p,p') = cj 21 (p,p') = l--(l + p-p') 2 , (A9) 



w 22 (p,p') - -(1 + p-p') 2 . 



(A10) 



APPENDIX B: ANGULAR INTEGRATIONS IN 
CSL PHASE 



In the calculation of the A-rate in the CSL phase, there 
are four different types of angular integrations over the 
phase space of quark momenta. Thus, the results for the 



functions, formally defined by Eq. (|43p in the 
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main text, have the following general structure: 



Aim 


= A 


+ K X 


+ K 2 


+ A" 3 , 


(Bla) 
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-K x 


+ K 2 


-ATs, 


(Bib) 


■P1121 


= K 


+ K 1 


-K 2 


- A3 , 


(Blc) 
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= K 


-K x 


-K 2 


+ K 3 , 


(Bid) 


Ami 


= K 


+ K\ 


-K 2 


-K 3 , 


(Blc) 


Fyivi 


= K 


— K\ 


-K 2 


+ K 3 , 


(Blf) 


F\ 22 \ 


= K 2 


+ A 3 


5 




(Big) 


F\ 222 


= K 2 


-K 3 


5 




(Blh) 



A2111 = K -K 1 +K 2 -K 3 , (Bli) 

A2112 = K± + K 3 , (Blj) 

A2121 = K - K x - K 2 + K 3 , (Blk) 

A2122 = A'i - K 3 , (Bll) 

A2211 = K - Ax - K 2 + K 3 , (Blm) 

A2212 = A'i - , (Bin) 

A2221 = K 2 -K 3 , (Bio) 

A2222 = A" 3 , (Blp) 



where the four types of angular integrals are given by: 



A" = J dQi J dfl 2 J dQ 3 J d£h(l - Pi ' P2)(l - P3 • p 4 )5(pi + p 2 - p 3 - P 4 ), (B2) 

Ai = J J dfii / dft 2 /" dft 3 y dQ 4 (l - Pi ■ p a )(l - Ps ■ p 4 )(l + p 4 ■ Pi) 2 £(pi + P2 p 3 - P 4 ), (B3) 

K 2 = jJ d£li y dtt 2 J dil 3 J dfi 4 (l - Pi • p 2 )(l - p 3 • P 4 )(l + Pa • P2) 2 5(pi + p 2 p 3 - P 4 ), (B4) 

K 3 = ^JdQxj dCl 2 y dn 3 y df2 4 (l - Pl • p 2 )(l - p 3 • p 4 )(l + p 4 • pi) 2 (l + p 3 • p 2 ) 2 5(pi + p 2 - p 3 - p 4 ) . (B5) 



The result for Kn was obtained in Rcf. 123. It reads 



4?r 3 

K = 2 2 2 2 A (pi2,Pl2,P3 4 ,-P3 4 ) , (B6) 
PlP2 P3P4 



where p {j = \p t -pj\, Pi 3 = p t +Pj, and 



L {a, b, c, d) = 0(c - a)Q(d - 6)6(6 - c) J (c, 6, 6, d) + 9(o - c)0(d - 6) Jo(o, 6, 6, d) 

+ 0(a - c)0(6 - d)0(d - a) J (a, d, b, d) + 0(c - a) 0(6 - d) J (c, d, 6, d) , (B7) 

which is given in terms of 

J (a,b,c,d) = ( dP(c 2 - P 2 )(d 2 - P 2 )=c 2 d 2 {b-a)--{c 2 + d 2 )(b 3 -a 3 ) + -(b 5 - a 5 ). (B8) 

Ja 3 5 

To leading order in powers of large n, this result simplifies to 

i (0,2/i,0,2 At ) = -|-. (B9) 
15 

By making use of Eq. (|B6|) . therefore, we obtain 



4tt 3 2 10 tt 3 
/if ~_i (0 J 2/i,0 > 2/i) = 1 ^ g -. (BIO) 

Using the same approach, in the following subsections wc calculate the results for K x , K 2 , and K 3 . 
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1. Calculation of K\ 

Here we calculate the angular integral K\ . Following the approach of Ref. [25| , we obtain 

Ki = i J dfli J dQ 2 J dfl 3 J dfl^l - pi • p 2 )(l - p 3 • p 4 )(l + P4 • Pi) 2 ^(Pi + P2 P3 - P4) 



4pi 



4 P i 



dfl x / dfl 3 / dQ 4 



1-— (Pl-P-Pl) 

P2 



#1 



d#i sin 0i 



(1 - p 3 • P4 )(l + P4 • Pl) 2 S{ P 2 - |P - Pl|) 

(1-P3-P4)5(P2-|P-P1| 



1 (P cos 6*i — pi 

P2 



[1 + cos 9i cos 4 + sin 9\ sin # 4 cos(<^i 



(BH) 



where P = P3 + P4 and P — |P|. In order to integrate over 9\, we choose the coordinate system so that the z-axis is 
along the vector P. After making use of the (5-function, we easily integrate over 9\ and arrive at the following result: 



K\ = / dfl 3 / dfl 



p2 p2 
^12 r 



(i - P3 • p 4 )e(p 12 - p)o(p - Pl2 ) 



sppIpI 

I d(j)i[l + cos9\ cos #4 + shi9\ sin(?4cos(0i — ^4)] 
Jo 

p2 p2 
M2 r 



7T dfl 3 dfl 



SPpfpi 



0(Pi2 - P)S{P - p 12 ){l - p 3 • P4) 2(1 + cos (9* cos <9 4 r + (sine* sin* 



(B12) 



where cos 6** = (P 2 + p 2 — p 2 )/2Ppi and cos04 = (p>4 • P) = (P 2 +p 2 — p 2 )/2Pp 4 . It may be appropriate to emphasize 
that the result of the last integration was presented in a form that independent of a specific choice of the coordinate 
system. As can be easily checked, the integrand in Eq. (|B12|) depends only on the relative angle 6*34 between the 
vectors P3 and P4 (or equivalently only on the variable P). Therefore, while integrating over we could fix the 
orientation of p>3 arbitrarily. It is convenient to choose P3 as the z-axis and perform the integration over The 
result is independent of the angular coordinates in Q3. Thus, the remaining integration over O3 gives an extra factor 
Air. In the end, we arrive at 



Ki = 



^pIpIpIpI 



p 3 1 



g(P)dP, 



(B13) 



1>3 1 



where 



g(P) = {P 2 2 -P 2 ){P 2 A -P 2 )Q(P 12 -P)Q{P~p 12 ) 



2 1 



(p 2 +pi-p 2 2 )(p 2 +pj-pi r 2 
AP 2 PlPi 



{p 2 +pI-p 2 2 ? 

AP 2 p\ 



{p 2 +pI-p I? 
ap 2 pI 



(B14) 



Note that we changed the integration variable from 6*34 to P = \J p\ + p 2 + 2p 3 p 4 cos $34. 

The final result for K\ can be conveniently given in the same form as Kq in the previous section, i.e. 



ip{p 2 p{pi 



Ll(pi 2 ,Pl 2 ,P34,P34) , 



where, by definition, 

Li(a, b, c, d) 



Q(c - a)e(d - 6)0(6 - c) J^c, 6, 6, d) + 9(a - c)Q(d - 6) Ji(a, b, 6, d) 
Q(a - c)Q(b - d)e{d - a) J^a, d, 6, d) + Q(c - a)9(6 - d) Ji(c, d, b, d) , 



(B15) 



(B16) 



and 



Ji(a,b,c,d) = 



dP{c 2 - P 2 ){d 2 - P 2 ) 
(P 2 +p 2 -p 2 ) 2 



2 1 



+ 1- 



AP 2 p\ 



| (P 2 +pI~ p 2 2 )(p 2 +pI-p\) 

AP 2 p lP4 

(p 2 +pI~pI) 2 ' 

AP 2 p\ 



(B17) 
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To leading order in powers of large [i, this result reduces to 



Li(0,2/i,0,2/i) = 



35 



This, in turn, gives 



7T 3 2 10 7T 3 



(B18) 



(B19) 



2. Calculation of A'2 

As is easy to see, the expression for K2 can be obtained from K\ by the following exchange of variables: p\ <-> p 2 
and j)3 <-> P4. Thus, the result reads 



Lz(pi2, Pl2, P34, -P34) , 



where 



L 2 (a, 6, c, d) = 0(c - a)0(d - 5)9(6 - c) J 2 (c, 6, 6, d) + 0(o - c)0(d - 6) J 2 (a, 6, 6, d) 
+ 0(a - c)0(6 - d)0(d - a) J 2 (a, d, 6, d) + 0(c - a) 0(6 - d) J 2 (c, d, 6, d) , 



and 



J 2 (a,b,c,d) = I dP(c 2 - P 2 ){d 2 - P 2 ) 



( (p'+p^-pIKp'+pI-pI) 

+ 4P 2 P2P3 



1 {p' 2 +pj-p{) 2 \ ( 1 (P 2 +Pi-Pl) 2 



4P 2 pl J \ \P 2 p\ 



We also find that Ki is identical to K\ to leading order in powers of large fx, i.e., 

2 10 7T 3 



35/i 3 



(B20) 



(B21) 



(B22) 



(B23) 



3. Calculation of K3 



Now we calculate the angular integral K3. We start by using the same approach as in the calculation of K\, 



l^f 1 



1-— (pi-P-.Pi) 

P2 



(1 -P3 " p 4 )(l +P4 ' Pi} 



P2P3 



(p 3 • P - Pi ■ P + Pi p 4 ) 



*(P2-|P- P 1| 



(B24) 



To calculate the integral over f2i, we fix the coordinate system so that the z-axis coincides with the direction of P. 
After integration, we obtain 



df2 4 - 



P\\ - P ' 



e(Pi 2 -p)6(p-pi 2 )(i-P3-P4) 



j d ^j—,2Pp 2 p 2 

/ d(j>i [1 + cos 8\ cos #4 + sin 9\ sin #4 cos(^i 
Jo 



, , P2-P1+P3-P4 , PiPi , a * a , ■ a* ■ a (a aw 
1 H — H (cost^i cos 04 + smfc>i sin 6*4 cos(0i — 04)1 

2p2P3 P2P3 



(B25) 
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where cos 6** = (P 2 + p\ — p?,)/2Ppi and cos#4 = (p4 • P) = (P 2 + p 2 — p 2 )/2Pp4. In the derivation, we also used the 
following relation: 



P3 P 

By performing the integration over <f>i , we derive 

[Ph-P 2 ){Ph-P 2 ' 



p 2 +pI-pI 



(B26) 



K 3 = 7T / dfl 3 / dn 



{Pi -Pi,f 



±PiPi {pL - Pi 3 ) i 



imPpipipM e ( P 12 -P)e ( P-, 12 ) 

/ {P 2 +pl-p 2 2 ){P 2 +pl~pi) \ 2 
+ 4P 2 PlPi 1 

(p 2 +p 2 1 -p 2 2 )(p 2 +pI-pI) 



{p 2 +pI-pI? 
\p 2 p\ 



i - 



(p 2 + pI-pD 2 

4P 2 p 2 



iP 2 p lPi 

2 i „2 „2Wp2 i „2 ^,2 \ \ 2 



2 ^ , 



2 2 



-24 



(P 2 +^-pD(P 2 +^-pi) 

4P 2 piP4 



( (P 2 +p 2 -p 2 ) 2 \ ( (p 2 +pI-p 2 3 ) 2 \ 

\ ±P 2 Pl A 4P2p 2 ^ 

32 1 ^,2 ^2\2N 2 



3 1- 



(p 2 +pj- P 2 2 y 

4F 2 p 2 



1 - 



(p 2 + P i-p 2 3 y 

±p 2 pI 



4P 2 p lPi 



AP 2 p\ 



1 , (p 2 +p 2 -^)(p 2 +p!-p1) V/ 1 (p 2 +p 2 -pI) 2 A A (p 2 +p!~p 2 ) 2 



4P 2 p 2 



(B27) 



Finally, in order to integrating over we use the coordinate system with the z-axis along p3. Then, we get 



Z2p\p\p%p\ 
where 

HP) = (P 1 2 2 -P 2 )(P 3 2 4 -P 2 )e(P 12 -P)6(P-p 12 ) 

?2 D 2 \ 2 



P34 



h{P)dP, 



(B28) 



(A 2 4-P 2 2 3 )' 



4P1P4 (Pf4 - P 2 2 3 ) 1 



2 ( 1 (p 2 +pI~p 2 2){p 2 +pI-pD 

4P 2 P1P4 

(P 2 +p 2 1 -p 2 2 )(P 2 +pl-p 2 3 ) 



1 - 



(P 2 +P 2 -P 2 ) 2 

AP 2 P 2 



1 - 



(p 2 +pI-pI) 2 
ap 2 p\ 



4P 2 p lP4 



2 ( 1 , {p 2 +pI-pI){p 2 +pI-pI) 

4P 2 p lP4 



AP 2 p 2 



1 {p 2 +p\-pI) 2 \ ( 1 (p 2 +pI-pD 2 



AP 2 p 2 



2 2 
PiPi 



1 , (P 2 +P 2 -p|)(^ 2 +pl-pI) y , 2 ( 1 ( p2 +p\-P2) 2 ) 2 d (p 2 +pI-pD 2 



4P 2 p lP4 

j2 , _ „2\(r>2 , m 2_ m 2 



AP 2 p 2 



+ 94 (, , (^ 2 +p 2 -p 2 .)op 2 +p 2 -p 2 ) v / _ (p^^in / _ {p 2 +pi-pd- 

\ 4P 2 pipa ) \ ±P 2 p\ J \ 4P 2 p 2 

The result can be given in the same form as the other integrals in the previous subsections, i.e., 



(B29) 



^p\p\p\p\ 



£3(^12, Pl2,P34,P34) , 



(B30) 



where 



L 3 {a, 6, c, d) = 9(c - a)@(d - 6)9(6 - c) J 3 (c, 6, 6, d) + 6(a - c)6(d - 6) J 3 (a, 6, 6, d) 
+ 6(a - c)6(6 - d)6(d - a) J 3 (a, d, 6, d) + 6(c - a)9(6 - d) J 3 (c, d, 6, d) , 



(B31) 



16 



and 

J 3 (a,b,c, d) 



[ dP (c 2 - P 2 )(d 2 P 2 ){(P? 4 P| 3 ) 2 



2 1 



(P 2 + p \-pI){p 2 +pI-pI) 



1 - 



(p 2 +p\-pI? 



(i? 4 - P&) ( 1 



4P 2 pip 4 

, (p 2 +pI-pI){p 2 +pI-pI) 



4P 2 p 2 



1 - 



(p 2 +pI-pI) 2 

AP 2 p 2 



4P 2 p lP4 



2 ( 1 | {p 2 +pI-pI){p 2 +pI-pI) \ 

\ 4P 2 PlP4 J 



1 {p 2 +p\-pD 2 \ ( 1 _{p 2 +pI-pI? 



2 2 

pipa 



1 ! (P 2 +p 2 -pl)(P 2 +pI-pj) 
4P 2 pip 4 

94 A , (p 2 +p 2 -^)(p 2 +p!-pI) \ 2 



3 1- 



4P 2 p 2 
(P 2 +p 2 - 

4P 2 p 2 



1 - 



4P 2 p 2 

(p 2 + P i-Pir 2 

AP 2 p 2 



AP 2 p 2 



1 (P 2 +Pi-Pi) 2 \ ( l (p 2 +pI-pI) 2 



AP 2 p 2 



(B32) 



r 



To leading order in powers of large fj,, the result reduces 
to 



L 3 (0,2 M ,0,2 M ) 



29 x 2 2 V 
45045 



(B33) 



Then, by making use of the relation in Eq. (|B30|) . we 
derive 



4 F 



2/' 



12 



L3(0,2/x,0,2/x) 



29 x 2 15 tt 3 
45045/i 3 ' 



(B34) 



to leading order in inverse powers of ji 



By making use of the leading order results for Ki 
(i = 1, 2, 3) obtained in the previous subsections, here we 
write down the explicit results for the functions Fr 1 r 2 r 3 r l 
in the same approximation: 



1301 x 2 10 tt 3 
9009m 3 



(B35a) 



-F1221 
P1112 
P1122 
P1222 
P2222 = 

Note that 



P2112 
P1121 
P212 



443 x 2 1, V d 
9009/i 3 



(B35b) 



415 x 2 10 tt 3 , 
P1211 = Fam = - 5 , (B35c) 



P2121 — P2211 — 



P2122 — P2212 — P2221 
29 x 2 15 tt 3 



9009/x 3 

1357 x 2 10 tt 3 
45045^3 

359 x 2 10 tt 3 
45045^ 3 ' 



45045/i 3 



,(B35d) 
(B35e) 
(B35f) 



^ ri7 



ri,r 2 ,r 3 ,r 4 



3 x 2 1( V 
5 M 3 ' 



(B36) 
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